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Introduction 


The  clinical  course  of  breast  cancer  is  highly  variable  between  patients,  particularly  when 
it  is  discovered  at  its  earliest  stage.  For  example,  of  patients  with  newly  diagnosed  Stage  I  breast 
cancer,  approximately  20%  will  subsequently  develop  metastatic  disease  within  five  years  while 
80%  will  remain  disease  free  after  resection  alone  (1).  Because  of  the  controversial  cost  and 
morbidity  to  total  benefit  ratio  of  treating  all  Stage  I  breast  cancer  patients  with  aggressive 
chemotherapy,  an  enormous  emphasis  has  been  placed  on  the  development  of  biological  markers 
to  better  assist  in  treatment  planning  and  prognosticating  the  clinical  course  of  breast  cancer 
patients  (2).  In  addition  to  traditional  analyses  of  tumor  size,  lymph  node  involvement,  hormone 
receptor  function,  and  nuclear  grade,  a  plethora  of  molecular  markers  have  been  proposed  as 
prognostic  tools  (3-7).  The  biological  rationale  supporting  the  clinical  utility  of  many  of  these 
markers  is  that  the  evolution  of  breast  cancer  is  accompanied  by  multiple  changes  in  gene 
expression  and  that  quantification  of  these  changes  can  provide  predictive  information  about  the 
clinical  course  of  individual  tumors.  However,  while  analysis  of  individual  gene  expression  has 
been  shown  to  have  predictive  value  within  single  studies,  there  is  no  clear  consensus  to  suggest 
that  analysis  of  any  one  gene  or  gene  product  has  independent  predictive  power  (8). 

As  our  understanding  of  both  the  human  genome  and  the  molecular  genetic  basis  of 
cancer  has  expanded  over  the  past  five  years,  it  is  now  believed  that  a  broader  examination  of 
global  patterns  of  gene  expression  may  be  a  more  sensitive  and  specific  method  for 
characterizing  tumors  at  the  molecular  level.  A  number  of  technical  advances  have  made  such 
an  approach  possible.  Catalogs  of  thousands  of  expressed  genes  in  the  form  of  Expressed 
Sequence  Tags  (ESTs)  are  now  available  for  analysis  (9-1 1).  Many  of  these  catalogs  are  derived 
from  tumor  tissue  sources,  including  breast  cancer.  The  availability  of  this  sequence  information 
enables  investigators  to  identify  potential  genes  that  contribute  to  the  pathogenesis  of  cancer  and 
that  may  ultimately  serve  as  prognostic  markers.  A  second  technical  advance  is  the  ability  to 
simultaneously  and  quantitatively  measure  expression  of  multiple  genes  using  cDNA  arrays  (12- 
14).  This  experimental  approach  allows  for  the  analysis  of  global  patterns  of  gene  expression 
rather  than  traditional  methods  that  examine  expression  ‘one  gene  at  a  time’  (15). 

While  both  of  these  strategies  represent  important  technical  advances,  neither  has  been 
specifically  directed  toward  addressing  clinical  problems.  For  example,  EST  catalogs  from 
tumor  sources  (such  as  breast  cancer)  contain  sequences  that  are  neither  tumor-specific  nor  cell- 
specific.  Therefore,  identification  of  clinically  relevant  gene  sequences  will  require  considerable 
filtering  of  this  sequence  information.  Furthermore,  the  EST  sequences  obtained  from  tumor 
specimens  contain  no  accompanying  information  with  regard  to  the  tumor  or  patient  from  which 
they  were  obtained.  For  this  reason,  it  will  be  impossible  to  use  this  sequence  information  to 
correlate  patterns  of  gene  expression  with  any  clinically  or  pathologically  significant  parameters. 
While  still  in  its  infancy,  the  development  of  cDNA  arrays  for  the  multiplex  analysis  of  gene 
expression  has  emphasized  the  placement  of  large  numbers  of  randomly  chosen  cDNAs  in  an 
array  to  build  a  comprehensive  screen  of  gene  expression.  Without  a  rationale-based  approach 
for  selecting  gene  sequences  to  place  in  the  probe  array,  a  very  large  array  will  be  needed  and  the 
majority  of  patterns  of  expression  detected  will  be  irrelevant  to  individual  tumor  phenotypes. 

Our  interest  is  in  identifying  quantitative  changes  in  gene  expression  that  correlate  with 
significant  clinical  parameters  of  human  breast  cancer.  By  identifying  a  finite  set  of  these 
changes,  we  believe  that  the  patterns  of  gene  expression  (‘expression  fingerprints’)  of  individual 
tumors  will  be  a  sensitive  and  specific  diagnostic  tools  for  assisting  in  prognosis  and  therapeutic 
management  of  breast  cancer  patients.  We  have  previously  demonstrated  that,  using  several 
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modifications  of  the  differential  display  polymerase  chain  reaction  technique  (DDPCR),  we  are 
able  to  create  a  high  throughput  system  that  can  rapidly  identify  Differentially  Expressed 
Sequence  Tags  (DESTs)  from  primary  human  breast  tumors  (16).  These  DESTs  are  cDNA 
fragments  corresponding  to  mRNAs  that  are  differentially  expressed  in  breast  tumors  as 
compared  to  patient-matched  normal  breast  tissue.  A  catalog  of  these  sequences  would  provide  a 
potential  list  of  new  molecular  markers  for  phenotyping  breast  tumors,  particularly  if  they  are 
correlated  with  clinical  information  about  the  tumor  and  patient  from  which  they  were  derived 
(16).  In  fact,  we  have  already  demonstrated  that  one  DEST  obtained  in  this  manner  is  a 
potentially  useful  new  molecular  marker  for  breast  cancer  (17). 

In  this  report,  we  describe  our  progress  in  adopting  our  preliminary  work  toward  the 
creation  of  a  breast  tumor  DEST  library,  similar  in  scope  to  EST  libraries  that  are  now  in 
existence,  but  with  added  clinical  relevance  to  human  breast  cancer.  Thus  far,  we  have 
optimized  protocols  using  five  different  human  breast  tumor  cell  lines  to  isolate  150  DEST 
fragments.  Having  demonstrated  the  feasibility  of  this  approach  using  human  breast  tumor  cell 
lines,  we  are  now  in  the  process  of  analyzing  a  small  set  of  well  defined,  microdissected,  primary 
human  breast  tumors.  We  believe  the  results  of  this  work  will  ultimately  define  a  set  of  gene 
sequences  whose  pattern  of  expression  could  collectively  define  a  breast  tumor  with  regard  to  its 
metastatic  potential,  therapeutic  vulnerability,  and  prognostic  significance  to  the  patient. 

Experimental  Methods 

The  differential  display  polymerase  chain  reaction  (DDPCR)  and  RNA  fingerprinting 
have  been  used  extensively  to  isolate  genes  that  are  differentially  expressed  between  different 
cell  and  tissue  types  (18, 19).  In  the  DDPCR  technique,  messenger  RNA  (mRNA)  is  isolated 
from  multiple  cell  sources  and  reverse  transcribed  to  produce  a  set  of  cell-specific 
complementary  DNAs  (cDNAs).  These  cDNAs  are  then  used  as  template  for  a  polymerase 
chain  reaction  (PCR)  at  low  stringency  annealing  conditions.  Polymerization  is  primed  with  one 
oligonucleotide  of  an  arbitrarily  defined  sequence,  ten  nucleotides  in  length  (lOmer),  and  a 
TuMN  oligonucleotide  (M=A,G,C;  N=A,C,G,T).  In  the  presence  of  radiolabeled  nucleotides, 
the  resulting  randomly  primed  polymerization  products  theoretically  represent  a  subset  of 
expressed  mRNAs  from  each  cell  population.  Each  PCR  reaction  is  run  side  by  side  on  an 
acrylamide  sequencing  gel  and  products  unique  to  one  cell  type  are  identified,  isolated,  and 
subcloned  for  further  study.  We  have  made  a  number  of  innovative  modifications  to  this 
technique  to  increase  the  yield  of  truly  differentially  expressed  sequences  that  are  isolated,  and  to 
allow  for  the  rapid,  high-throughput  analysis  of  these  sequences  without  the  time-  and  labor- 
intensive  process  of  subcloning.  Using  this  improved  method,  we  have  identified  differences  in 
gene  transcription  among  five  different  human  breast  tumor  cell  lines  (Table  I).  Figure  1 
outlines  this  procedure. 

Indicated  cell  lines  were  obtained  from  the  American  Type  Culture  Collection  and  grown 
in  media  formulations  as  recommended  by  the  supplier.  Approximately  lxlO6  cells  from  each 
cell  line  were  isolated  in  duplicate.  From  each  pair  of  cell  samples,  5  pg  of  total  RNA  was 
purified,  treated  with  RNase-free  DNase,  reverse  transcribed,  degraded  with  RNaseH  and 
DNase-free  RNase,  and  purified  by  spin  chromatography.  After  considerable  experimentation 
with  different  cDNA  synthesis  and  purification  schemes,  we  found  that  these  steps  of  first  strand 
cDNA  preparation  were  absolutely  necessary  to  produce  reliable  and  reproducible  DDPCR 
displays. 
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Synthesized  cDNA  pools  from  each  cell  line  were  subjected  to  the  DDPCR  technique  as 
described  above  and  previously  (18, 19)  except  for  several  notable  changes.  First,  we 
investigated  the  use  of  several  different  32P  radiolabeled  Tmer'  oligonucleotides  to  generate 
labeled  PCR  products.  We  compared  a  set  of  four  degenerate  TnVN  oligonucleotides  (where  V 
is  a  mix  of  A,  G,  and  C  bases  and  N  is  either  A,  C,  G,  or  T),  a  set  of  three  specific  T i  iN 
oligonucleotides,  and  a  set  of  six  (of  a  total  possible  of  12)  specific  TnNN  oligonucleotides. 
After  optimizing  annealing  temperatures  and  cycling  conditions  for  each  primer  set,  we 
concluded  that  the  set  of  four  degenerate  TnVN  oligonucleotides  produced  the  most  specific  and 
reliable  display  pattern.  Each  of  the  TnVN  oligonucleotides  also  contained  a  10  nucleotide  5' 
'anchor'  sequence  corresponding  to  the  T7  RNA  polymerase  promoter  for  the  purposes  of 
secondary  amplification  of  isolated  PCR  fragments.  End-labeling  oligonucleotides  with  T4 
polynucleotide  kinase  (as  opposed  to  incorporation  of  32P-labeled  dCTP)  produced  superior 
results  in  terms  of  pattern  reproducibility  and  theoretically  ensured  that  all  labeled  (and  hence 
visualized)  PCR  products  were  the  result  of  priming  from  mRNA  polyA  tails.  The  second 
(lOmer)  primer  employed  in  the  PCR  reaction  was  an  arbitrary,  ten  nucleotide  sequence  with  an 
additional  ten  nucleotide  5'  'anchor'  sequence  corresponding  to  the  Ml  3  origin  of  replication 
sequence.  A  number  of  different  arbitrary  sequences  were  incorporated  in  the  'lOmer'  primer,  all 
of  which  were  devoid  of  palindromic  sequences  and  maintained  a  G/C  content  of  60%. 

Although  the  pattern  of  amplified  fragments  differed  using  each  'lOmer'  primer  (see  figure  2), 
the  overall  performance  of  the  reactions  were  not  significantly  affected  by  primer  sequence 
composition.  After  considerable  trial,  a  set  of  four  degenerate  'Tmer'  and  seven  'lOmer'  primers 
were  established  for  all  future  work  (Table  II).  Finally,  conditions  were  optimized  such  that 
reactions  could  be  reliably  performed  in  a  10  pi  volume  and  loaded  directly  on  a  4%  denaturing 
acrylamide  sequencing  gel  for  preparative  analysis. 

For  one  given  experiment,  duplicate  first  strand  cDNA  reactions  were  performed  with 
each  human  breast  tumor  cell  line  as  described  above.  A  single  display  experiment  involved 
amplification  and  electrophoresis  of  10  duplicated  pairs  of  cDNA  samples  using  a  single  'Tmer' 
primer  in  combination  with  three  different  'lOmer'  oligonucleotides  (10x1x3=30  total  samples). 
Each  experiment  was  repeated  in  triplicate  to  verify  the  reproducibility  of  the  display  pattern.  A 
representative  display  is  shown  in  figure  2.  Experiments  were  then  repeated  for  each  of  the 
seven  'lOmer'  primers  permuted  with  each  of  the  four  'Tmer'  sequences  described  above.  Unlike 
other  studies  described  to  date,  potentially  relevant  DEST  fragments  were  selected  for 
reamplification,  sequencing,  and  entry  into  the  DEST  database  using  a  set  of  strict  criteria. 

These  criteria  were:  i)  Fragment  absent  (loss  of  expression)  in  duplicate  cDNA  reactions  from  at 
least  two  different  cell  lines  or  ii)  Fragment  present  (gain  of  expression)  in  duplicate  cDNA 
reactions  from  at  least  one  different  cell  line.  In  addition,  fragments  chosen  for  further  analysis 
could  not  be  present  in  the  same  cell  line  cDNA  reactions  that  employed  two  different  'lOmer' 
primers.  This  criteria  ensured  that  the  amplified  band  was  synthesized  by  the  priming  of  both 
the  specific  'lOmer'  primer  and  specific  'Tmer',  rather  than  the  'Tmer'  alone.  Bands  resulting 
from  'Tmer'-only  priming  events  are  not  candidates  for  direct  PCR  sequencing  since  their  ends 
are  symmetrical. 

Candidate  bands  were  excised  from  their  acrylamide  sequencing  gel,  eluted,  and 
reamplified  in  a  standard  polymerase  chain  reaction  using  the  original  'lOmer'  and  'Tmer' 
oligonucleotide  primers  employed  for  DDPCR.  We  examined  several  other  strategies  that 
employed  nested  primer  approaches  to  affect  secondary  amplification,  but  found  that  this  was 
unnecessary.  Reamplified  bands  were  quantitated  and  qualitatively  assessed  by  agarose  gel 
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electrophoresis  (figure  3).  Approximately  20  ng  of  each  secondarily  amplified  PCR  product  was 
subjected  to  cycle  sequencing  using  a  32P  end-labeled  oligonucleotide  corresponding  to  the 
'lOmer'  used  for  the  original  DDPCR  reaction.  We  again  examined  the  use  of  nested  sequencing 
primers  to  enhance  sequencing  read,  but  found  that  this  strategy  did  not  improve  results. 

Isolated  sequence  information  was  compared  to  the  NCBI  Genbank  database  (20,  21)  to 
determine  whether  a  DEST  fragment  represented  a  previously  characterized  gene,  and  was  then 
entered  into  a  computer  database  designed  in  Microsoft  Access  format  (figure  4). 

Results  and  Discussion 

Using  the  methodology  described  above,  we  have  attempted  to  isolate  gene  sequences 
that  are  differentially  expressed  among  five  different  human  breast  tumor  cell  lines.  Figure  2 
demonstrates  a  sample  result  from  the  numerous  differential  display  experiments  that  we  have 
performed  over  the  past  year.  We  have  identified  148  PCR  fragments  that  are  reproducibly  and 
differentially  present  among  the  cell  lines  tested.  Figure  3  demonstrates  a  sample  result  from 
isolation  and  secondary  amplification  of  a  small  set  of  these  fragments.  Overall,  approximately 
90%  of  bands  identified  and  isolated  resulted  in  secondary  amplification  products.  We  have 
used  a  direct  cycle  sequencing  approach  to  sequence  and  identify  the  first  set  of  these 
secondarily  amplified  products.  Thus  far,  approximately  70%  of  products  are  amenable  to  direct 
sequencing.  Figure  5  demonstrates  sample  sequencing  reads  from  several  of  these  PCR 
products.  Although  the  overall  quality  of  the  sequence  is  poor,  we  can  routinely  read  50-100 
nucleotides  of  sequence  with  ~95%  accuracy.  This  is  sufficient  sequence  information  to 
uniquely  identify  the  expressed  gene  and  to  design  tertiary,  gene-specific  primers  for  future  use. 
Each  sequenced  DEST  has  been  entered  into  a  database  that,  upon  further  maturation,  will  be 
made  available  to  the  public  realm  via  internet  access.  Figure  4  demonstrates  a  sample  data 
entry  form  in  this  database.  Although  we  are  still  only  beginning  to  analyze  sequence  data, 

Table  III  list  several  differentially  expressed  genes  that  we  have  identified  thus  far.  The  fact 
that  at  least  three  of  these  sequences  have  already  been  identified  in  independent  experimental 
runs  suggests  that  the  sequences  that  we  are  isolating  are  not  artifacts  of  the  DDPCR  assay. 

Recommendations 

In  the  original  proposed  plan  of  research,  the  first  six  months  of  work  were  to  be  devoted 
to  creating  a  library  of  Differentially  Expressed  Sequence  Tags  (DESTs)  from  12  initial  cases  of 
human  breast  tumor/normal  pairs.  To  simply  the  system  in  which  we  intended  to  identify  and 
correct  technical  problems,  we  have  chosen  to  focus  on  human  tumor  cell  lines  that  can  provide  a 
more  abundant  and  controlled  source  of  RNA.  Furthermore,  despite  our  familiarity  with  the 
DDPCR  approach,  we  encountered  a  number  of  technical  difficulties  that  required  re¬ 
optimization  of  protocols.  As  such,  Task  1  has  been  extended  through  month  12  of  the  proposal. 
Also  in  the  original  plan,  months  7-12  were  to  be  devoted  to  establishing  an  initial  cDNA  array 
of  DEST  sequences  and  developing  protocols  to  quantitate  gene  expression  using  radiolabeled 
cDNAs  from  2-5  sets  of  human  tumor/normal  tissue  pairs.  Because  of  the  delay  in  establish  a  set 
of  candidate  probe  sequences,  this  task  has  not  been  initiated. 

In  view  of  our  current  progress  and  stated  aims  in  the  original  proposal,  we  propose  the 
revised  time  line  for  this  study: 
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Task  1:  Months  13-14:  Complete  sequence  analysis  and  identification  of  148  sequences  derived 
from  human  breast  tumor  cell  lines.  Optimize  remainder  of  direct  cycle  sequencing  protocol  for 
secondarily  amplified  PCR  products. 

Task  3:  Months  15-24:  Prepare  duplicate  RNAs  and  first  strand  cDNAs  from  15  cases  of 
histopathologically  matched  breast  tumors  using  laser-mediated  microdissection.  Use  these 
cDNAs  to  isolate  at  least  200  different  DEST  sequences  with  40  combinations  of  primer  pairs. 
Enter  DEST  sequences,  clinical  data,  and  pathology  data  into  DEST  database.  The  objective  will 
be  to  use  the  now  established  experimental  approach  to  identify  molecular  variability  among 
otherwise  histologically  indistinguishable  breast  tumors  and  correlate  these  differences  to  other 
significant  clinical  variables. 

Task  4:  Months  15-20:  Synthesize  RNA  standards.  Create  preliminary  array  of  44  DEST 
sequences  identified  in  Task  1.  Test  hybridization  conditions  and  quantification  of  gene 
expression  using  radiolabeled  cDNAs  from  the  original  five  breast  cell  lines  used  in  Task  1. 
Identify  technical  problems  and  modify  protocol.  Evaluate  alternate  approaches  to  quantitate 
DEST  expression  levels  in  tumor  samples  such  as  multiplex  PCR  assays  using  semi-quantitative 
or  quantitative  (real-time)  assay  formats. 

Task  5:  Months  21-24:  Use  assay  format  established  in  Task  4  to  quantitate  DEST  gene 
expression  in  20  sets  of  primary  human  breast  tumors.  Enter  patterns  of  expression,  clinical  data, 
and  pathology  data  into  DEST  database. 


Conclusions: 

Although  progress  has  been  slower  than  anticipated  due  to  a  number  of  technical 
difficulties,  we  believe  that  the  current  study  may  still  result  in  a  much  larger  effort  directed  at 
defining  a  complete  set  of  diagnostic  expressed  sequence  tags.  For  example,  if  the  results  of 
these  experiments  continue  to  be  encouraging,  we  would  consider  using  even  newer  technical 
advances  to  display  a  more  complete  set  of  expressed  sequences  and  to  create  larger  arrays  of 
DEST  probes  (22).  We  would  also  expand  our  catalog  of  DEST  sequences  to  include  tumor 
specimens  with  accompanying  patient  data  describing  treatment  response,  disease-free  survival, 
and  overall  survival.  Just  as  we  are  attempting  to  correlate  DEST  expression  with  basic  clinical 
and  pathological  data  in  this  proposal,  we  could  eventually  correlate  expression  of  these 
sequences  with  this  more  clinically  pertinent  data.  Ultimately,  we  believe  that  this  approach  can 
be  used  to  create  a  nearly  complete  set  of  genetic  markers  whose  pattern  of  expression  would 
define  the  clinical  course  and  therapeutic  vulnerability  of  individual  patient  tumors.  Obviously, 
such  a  diagnostic  tool  would  have  tremendous  impact  on  the  diagnosis  and  treatment  of  breast 
cancer  patients. 
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FIGURE  1:  Scheme  for  Cataloging  DEST  Expression  Patterns  from  Human  Breast  Cancer 
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FIGURE  5:  Representative  Sequence  Analysis  of 
Reamplified  DEST  fragments.  Seven  secondarily 
amplified  DEST  fragments  were  subjected  to  direct  cycle 
sequencing  as  described  in  the  text.  Note  that  some 
fragments  yielded  very  good  sequence  ( 1 ,  6,  and  7), 
others  produced  poor  sequencing  reads  (3,  4),  and  others 
produced  no  readable  sequence  at  all  (2,  5).  Thus  far, 
approximately  70%  of  isolated  fragments  have  provided 
sufficient  sequencing  information  to  positively  identify 
fragments  as  either  previously  characterized  or  novel 
gene  products. 
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T47D 

MCF-7 

MDA-MB361 

MDA-MB415 

MDA-MB468 


Table  I.  Human  Breast  Carcinoma  Cell  Lines  Used  for  Generating  a  DEST  catalog 
for  Human  Breast  Cancer. 


Primer  Name 

Sequence 

T7TnVA 

5'-C  ACTAT  AGGGTTTTTTTTTTTV  A-3 '  | 

T7TnVC 

5'-CACTATAGGGTTTTTTTTTTTVC-3' 

T7TUVG 

5'-C  ACTAT AGGGTTTTTTTTTTTV  G-3' 

T7TnVT 

5'-C  ACTAT AGGGTTTTTTTTTTTVT-3' 

MF1 

5'-CAGTCACGACGACCGCTTGT-3' 

5'-C  AGT  C  ACG  AC  AGGTG  ACCGT -3’ 

MF3 

5'-CAGTCACGACTCCGGCTGAA-3' 

MF4 

5'-CAGTCACGACGCTGCGTGAT-3' 

MF5 

5'-CAGTCACGACGGTGCCTGAA-3' 

MF6 

5'-CAGTCACGACCTGCGGTGAT-3' 

5'-CAGTCACGACCGGCTGTGAA-3' 

Table  II.  Primer  Sequences  Used  for  Differential  Display  PCR. 


Cell  Line  Expression  Pattern 


DEST# 

Sequence  ID 

T47D 

MCF7 

MB415 

MB468 

MB361 

1 

ATP  Synthetase  p  Subunit 

- 

- 

+++ 

- 

2 

ATP  Synthetase  p  Subunit 

+ 

+ 

+++ 

+ 

+ 

3 

AA505271 

- 

- 

+++ 

- 

- 

4 

Ferratin 

- 

- 

+++ 

- 

- 

5 

AA992438 

+++ 

+++ 

+++ 

+++ 

6 

AA406625 

+++ 

- 

- 

- 

- 

7 

Ribosomal  Protein  L22 

+++ 

- 

- 

- 

- 

8 

ATP  Synthetase  p  Subunit 

- 

- 

+++ 

- 

- 

9 

Amplaxin 

+ 

+ 

+++ 

+ 

+ 

imu 

Amplaxin 

+ 

l 

+ 

+++ 

+ 

+ 

Table  III.  Representative  Selection  of  DEST  Sequences  Isolated  that  Correspond  to 
Previously  Known  Sequences.  Qualitative  expression  of  each  DEST  in  five  different 
human  breast  tumor  cell  lines  is  also  indicated. 
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